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PREFACE 


The  work  presented  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  at 
the  request  of  the  Ames  Research  Center  (ARC),  National  Aeronautics 
and  Space  Administration  (NASA),  under  Program  Element  65802F.  The 
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1.0  INTRODUCTION 


In  V/STOL  testing,  a partially  ventilated  test  section  (slotted  or 
perforated  walls)  offers  substantial  reductions  in  lift  and  blockage  inter- 
ference (Refs.  1 and  2).  However,  except  for  a few  cases,  only  mathe- 
matical models  typical  of  conventional  aircraft  (i.e.  sources  and  sinks 
or  doublets  for  blockage  interference  and  horizontally  trailing  horse- 
shoe vortices  for  upwash  interference)  have  been  employed  in  ventilated 
wall  analysis.  Lo  (Ref.  3)  simulated  a V/STOL  model  in  an  ideal  slotted 
tunnel  by  approximating  the  wake  originating  from  a rotor  by  means  of 
a skewed  cylinder  of  vortex  rings  of  constant  strength  lying  in  planes 
parallel  to  the  rotor  plane.  This  was  an  adaptation  of  the  mathematical 
model  suggested  by  Heyson  (Ref.  4)  and  Wright  (Ref.  5)  for  application 
in  closed  and/or  open  rectangular  tunnels. 

A recent  review  (Ref.  6)  of  the  state-of-the-art  of  low-speed  wind 
tunnels  for  V/STOL  testing  has  indicated  a need  for  better  mathematical 
simulation  of  the  V/STOL  model  with  a deflected  wake  in  theoretical  ap- 
proaches to  the  calculation  of  wind  tunnel  boundary  interference.  La- 
beled as  one  of  the  most  important  current  problems  was  the  curved 
wake  simulation  of  V/STOL  models  in  wall  interference  investigations. 
Heyson  (Ref.  7)  considered  the  effect  of  a uniform  strength  curved  wake 
in  test  sections  with  open  and/or  closed  boundaries  and  has  shown  that 
a correct  definition  of  the  wake  is  necessary  if  corrections  to  pitching 
moment  and  tail  forces  are  to  be  realistically  assessed.  Lo  (as  indi- 
cated in  Ref.  1)  performed  calculations  in  a tunnel  with  slotted  horizon- 
tal and  solid  vertical  walls  by  distributing  vortex  rings  on  an  empirically 
determined  curve  representing  the  centerline  of  a jet  in  crossflow. 

Binion  (Ref.  1)  compared  Lo's  theory  with  experimental  data  and  indi- 
cated the  lack  of  agreement  between  theory  and  experiment  may  be  as- 
sociated with  the  boundary  condition  for  a slotted  wall. 

Both  Heyson  and  Lo  assumed  that  the  jet  strength  along  the  curved 
wake  was  constant.  Kirkpatrick  (Ref.  8)  made  an  effort  to  include  the 
decay  of  the  jet  strength  along  the  jet  path  by  assuming  the  ring  vortex 
strength  to  be  proportional  to  the  cosine  of  the  wake  angle.  This  decay 
seemed  excessively  large  when  compared  with  the  rate  of  decay  of  a 
real  V/STOL  wake.  In  addition,  Kirkpatrick's  solutions  are  inaccurate 
because  limited  computer  capability  did  not  allow  the  use  of  a sufficient 
number  of  images  for  the  closed  wind  tunnel. 
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The  purpose  of  the  present  analysis  is  to  determine  the  wind  tunnel 
boundary  interference  on  a V/STOL  model  in  a rectangular  tunnel  with 
ventilated  horizontal  and  solid  vertical  walls.  The  mathematical  repre- 
sentation of  the  V/STOL  model  accounts  for  the  curvature  and  decay  of 
the  wake.  Solid  vertical  walls  are  selected  not  only  for  mathematical 
convenience,  but  also  because  Lo  and  Binion  (Ref.  9)  have  shown  that 
the  upwash  interference  is  insensitive  to  the  porosity  of  the  vertical 
walls. 


2.0  GENERAL  ANALYSIS 


The  flow  field  is  considered  three-dimensional  and  is  treated  as 
steady,  nonviscous,  incompressible,  and  irrotational  for  the  purpose 
of  determining  first-order  tunnel  wall  interference  corrections.  It  is 
also  assumed  that  the  velocity  perturbation  at  the  wall  caused  by  the 
body  and  the  perturbations  at  the  body  induced  by  the  wall  are  small 
compared  with  the  free- stream  velocity.  The  field  equation  of  an  in- 
viscid,  incompressible  fluid  in  terms  of  the  perturbation  velocity  po- 
tential ($)  is  the  well-known  Laplace  equation: 


2.1  BOUNDARY  CONDITIONS 

The  tunnel  geometry  to  be  considered  consists  of  solid  vertical  walls 
and  ventilated  horizontal  walls.  For  the  solid  vertical  walls,  the  boun- 
dary condition  of  no  flow  through  the  wall  is 

H = 0 atY  = ±b  (2) 


For  the  ventilated  horizontal  walls,  the  homogeneous  boundary  con- 
dition derived  by  Baldwin,  et  al.  (Ref.  10)  is  used: 


+ + v d2<E 

<3X  ~ R <?z  " ^ dxdz 


at  Z 


±h 


(3) 
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where  R is  a porosity  parameter  which  accounts  for  the  viscous  effects 
in  the  slots  and  must  be  determined  experimentally,  and  K is  a geomet- 
ric slot  parameter  derived  by  Chen  and  Mears  (Ref.  11)  for  a thin  wall 
as 


<«> 

where  d is  the  periodic  spacing  of  the  slots  and  a is  the  slot  width.  The 
open  area  ratio  a/d  is  assumed  to  be  small.  The  derivation  of  Eq.  (4) 
is  based  on  the  assumptions  that  both  the  velocity  normal  to  the  slot  and 
the  perturbation  from  the  mean  flow  are  small  compared  with  the  un- 
disturbed tunnel  velocity,  that  the  pressure  is  constant  across  the  open 
portions,  and  that  the  flow  angle  is  zero  on  the  solid  portions  of  the 
wall.  In  addition,  it  is  assumed  that  the  plenum  pressure  equals  the 
undisturbed  free- stream  pressure  and  that  walls  are  taken  to  be  straight 
and  parallel  with  constant  width  slots. 

Equation  (3)  contains  the  boundary  conditions  for  other  tunnel  walls 


as  limiting  forms, 

Solid  wall 

R -*  0 or  K "♦  oo 

4>z  = 0 

Free-jet 

K -*  0,  R o© 

<J>X  = 0 

Ideal  slotted  wall 

R - oo 

<t>x  + (1/R)4> 

Perforated  wall 

K ->  0 

<J>X  + (1/R)<D 

For  convenience  in  data  presentation,  a new  slot  parameter,  P = 

(1  + K/h)-1,  and  anew  porosity  parameter,  Q = (1  + 1/R)"*,  are  intro- 
duced so  that  the  interval  0 <L  P,  Q <.  1 represents  the  entire  range  0 <_K, 
1/R<®. 

The  boundary  condition  for  an  ideal  slotted  wall  tunnel  can  also  be 
written: 


$ + K«3>z  = 0 at  Z = ±h  (5) 

As  stated  in  the  derivation  of  the  general  equation  [Eq.  (4)],  a basic  as- 
sumption used  is  that  both  the  velocity  normal  to  the  slot  and  the  per- 
turbations from  the  mean  flow  are  small  compared  with  the  undisturbed 
tunnel  velocity.  However,  for  high-lift  V/STOL  models,  cross-flow 
velocities  in  the  test  section  may  become  appreciable;  thus  the  velocity 
near  the  slots  can  be  quite  high.  Therefore,  near  the  slotted  walls,  the 
quadratic  terms  of  the  Bernoulli  equation  may  no  longer  be  neglected  as 
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in  the  derivation  of  Eq.  (5).  A derivation  of  a new  boundary  condition 
based  on  a heuristical  approximation  to  the  quadratic  terms  of  the 
Bernoulli  equation  for  a high-lift  V/STOL  model  in  an  ideal  slotted  tun- 
nel is  introduced  in  Appendix  A.  The  fundamental  assumption  of  the 
derivation  is  that  the  higher  cross-flow  velocities  near  a slotted  wall 
caused  by  a lift- augmented  model  may  be  approximated  by  a uniform, 
constant  cross  flow  related  to  the  increased  lift  induced  by  the  lift  aug- 
mentation device.  The  resulting  quasi-linear*  boundary  condition  for 
an  ideal  slotted  wall  takes  the  form  at  the  horizontal  walls 


dX  ~ R <9Z  “ Ke  dXdl 


_0  at  Z = ±h 


(6) 


where  the  "pseudo  porosity  parameter"  (1/Re)  is  defined  by 


= tan  «o  <7a) 

where  aQ  is  the  zero  lift  angle  of  attack  of  the  lift- augmented  V/STOL 
model.  The  "effective  slot  parameter"  (Ke)  is  twice  the  classic  slot 
parameter  defined  by  Eq.  (4).  The  new  quasi-linear  boundary  condition 
has  exactly  the  same  form  as  the  general  ventilated  wall  boundary  con- 
dition, Eq.  (3).  Also,  viscous  effects  in  the  slots  can  be  incorporated 
in  the  quasi-linear  boundary  condition  by  redefining  the  "pseudo  porosity 
parameter"  as 


l l 

— = — tan  a, 
Re  R 


(7b) 


Hence,  a boundary  condition  of  the  general  form  of  Eq.  (3),  can  also  be 
used  for  high-lift  V/STOL  models. 

In  addition,  the  upstream  and  downstream  end  conditions 

= 0 at  X=— oc 


<?3>  „ v 

Jx  = 0 at  x = +oc 


(8) 


must  be  satisfied. 


*The  boundary  condition  is  called  quasi-linear  in  the  sense  that  the 
quadratic  cross-flow  velocity  is  approximated  in  a linear  fashion.  In 
the  mathematical  sense,  Eq.  (6)  is  a fully  linear  equation. 
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22  MATHEMATICAL  REPRESENTATION  OF  A V/STOL  MODEL 

The  augmented  lift  of  the  V/STOL  model  is  considered  herein  to  re- 
sult from  a rotor,  lifting  fan,  or  lifting  jet.  The  disturbance  potential 
of  a rotor,  lifting  fan,  or  lifting  jet,  is  represented  by  an  elliptic  vortex 
cylinder  sheet  following  the  path  of  the  wake.  Among  the  mathematical 
V/STOL  models  chosen  in  previous  investigations,  Heyson  (Ref.  12)  as- 
sumed that  the  wake  flows  in  a straight  line  intersecting  the  lower 
boundary  at  some  point  behind  the  model  and  then  flows  along  the  floor. 
Wright  (Ref.  5)  used  a similar  model  except  that  the  wake  was  assumed 
to  break  through  the  lower  boundary  and  to  have  no  further  influence  on 
the  flow  field.  Lo  (Ref.  3)  argued  from  physical  grounds  that  Wright's 
model  is  the  more  reasonable  of  the  two  models.  Furthermore,  while 
Heyson's  model  can  be  readily  applied  to  a closed  or  open  jet  where  an 
image  system  can  be  used  to  represent  the  influence  of  the  boundaries, 
in  the  present  generalized  analysis  where  an  image  system  does  not 
exist  for  a ventilated  wall,  severe  mathematical  complications  arise 
because  of  the  proximity  of  the  jet  wake  to  the  lower  boundary  in  the 
far  downstream.  Thus,  from  a physical  and  mathematical  standpoint, 

the  wake-boundary  interaction  assumption  of  Wright  will  be  used. 

/ 

A jet  in  cross  flow  which  is  used  herein  to  represent  the  model  wake 
has  two  essential  characteristics.  First,  the  path  of  the  curved  jet  de- 
pends primarily  on  the  initial  jet  to  free- stream  velocity  ratio  and  the 
initial  jet  deflection  angle.  Second,  the  strength  of  the  jet  decreases 
rapidly  as  the  jet  progresses  downstream  because  of  the  mixing  process 
between  the  jet  and  the  free  stream  and  the  viscous  action  within  the  jet. 
In  the  present  work,  the  jet  wake  is  treated  in  a linear  fashion  by  assum- 
ing that  the  tunnel  boundaries  have  no  influence  on  the  jet  trajectory. 
Hence,  the  jet  path  is  known  a priori  and  can  be  defined  by  existing  em- 
pirical descriptions  such  as  that  by  Margason  (Ref.  13).  Unfortunately, 
the  mathematical  description  of  the  jet  decay  has  not  been  described  in 
the  literature;  hence,  the  present  work  will  be  restricted  to  an  approxi- 
mate description  of  the  jet  decay.  For  purposes  of  analytical  integra- 
tion of  the  interference  potential  along  the  curved  jet  path  (see  Section 
2.4),  the  jet  path  is  represented  by  a series  of  straight  line  segments 
approximately  equivalent  to  the  curved  path,  as  shown  in  Fig.  1.  Fur- 
thermore, the  strength  of  the  jet  is  assumed  constant  over  each  straight 
line  segment  but  may  vary  from  segment  to  segment. 

The  model  is  assumed  to  be  mounted  at  the  center  of  the  test  section 
and  parallel  to  the  tunnel  centerline  since  Lo  (Ref.  3)  has  shown  the  in- 
terference factors  to  be  a weak  function  of  the  angle  of  attack.  The 
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Zn 


STRAIGHT  LINE 
APPROXIMATION 


Figure  1.  Mathematical  representation  of  V/STOL  model. 


vortex  sheet  is  made  up  of  a continuous  distribution  along  straight  line 
segments  of  circular  vortex  rings  lying  in  planes  parallel  to  the  model 
plane.  The  element  of  potential  (dfl^)  induced  at  a field  point  P by  a 
vortex  ring  of  strength  (dT/dS)dS  is 

d<Dni  = (dr'dS)dS //(cos0'  hrr2)dA  (9) 

where  r is  the  distance  from  an  element  of  surface  of  the  vortex  ring  to 
the  point  P,  and  S is  the  distance  along  the  jet.  By  using  the  small  per- 
turbation assumption,  the  variation  of  r over  the  surface  A enclosed  by 
the  vortex  ring  is  assumed  small  enough  such  that  the  angle  9 between 
the  vector  X and  the  vector  distance  is  constant.  Hence,  the  integra- 
tion of-  Eq.  (9)  is 

d<tm  = (dP'dS)dS(A  cos  O' in2)  (10) 

which  can  be  expressed  in  terms  of  X,  Y,  Z,  and  S for  the  nth  element 
by  using 

i = <*n-£2  + Y2  + (Zn-£)2 


10 


AEDC-TR  -74-51 


and 


cos  6 = A • r /Ar  = (Z  -£)/r 

n n n ^ n 

where 

£ = S sin  pB  , c = -S  cos£n 

It  follows  that 


zn  ~ C 

(d<Dm)  = (Adr/JS)„  — — as 

n 4jt  r_ 


Hence  the  disturbance  potential  of  the  nth  element  is 


(4>m)  = / (Adr/dS)n 


dS 


4 7T  rf 


where  the  location  of  the  nth  element  is  given  by 


(ID 


(12) 


n —1 


xn  = X - 2 Lj  Sin  0.  , n = 1 


,2, 


N 


Z + 2 Lj  cos  j3-  , n = 1,  2,  . . . , N 

j=o  ’ 1 


(13) 


\and  Ln  is  the  length  of  the  nth  segment.  The  geometry  of  the  length  of 
the  line  segments  is  shown  in  Fig.  2. 


2.3  METHOD  OF  SOLUTION 

By  assigning  the  free- stream  velocity  (U^)  and  the  tunnel  half- width 
(b)  as  the  characteristic  velocity  and  length,  the  mathematical  system 
can  be  normalized  by  defining 

x = X/b,  y = Y/b,  z = Z/b,  s = S/b,  £n  = Ln.'b 
6 = e/U^b,  y = r/L'^b,  F = K/h,  A = h/b 

Thus,  the  field  equation  in  normalized  coordinates  is 


i +il  + ii-V- 

• W dy2  5z  2T 


(14) 
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with  the  normalized  boundary  conditions 


dy 


0 at  y 


±1 


(15) 


and 


d<f>  , v p d2<t>  x 1 dd> 

dx  dxdz  R dr. 


0 at  z 


+A 


(16) 


<j>  = 0 at  x = — oo  ; dcfr-'dx  = 0 at  x = +* 


(17) 


The  linearity  of  the  field  equation  and  its  boundary  conditions  per- 
mits the  perturbation  potential  to  be  composed  of  two  parts  as 

<t>  - (18) 

where  <pm  is  the  potential  of  the  flow  about  the  model  in  free  air  and  <J>j 
is  the  interference  potential  induced  by  the  tunnel  boundaries. 
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If  0m  ,is  taken  to  be  a known  solution  of  the  Laplace  equation  which 
approximates  the  free- air  potential  at  points  far  from  the  model,  0^  can 
be  calculated  from  the  fact  that  the  sum  (0m  + 0j)  satisfies  the  known 
boundary  condition  at  the  wall.  Since  the  values  of  0m  are  used  only  at 
the  wall,  any  approximate  representation  (0m)  for  the  model  would  not 
appreciably  affect  the  calculation  of  0^. 

Since  the  model  potential  is  a series  of  vortex  rings  which  satisfy 
Laplace's  equation,  the  differential  equation  for  the  interference  poten- 
tial is 


V20;  = 0 (19) 

The  solution  of  Eq.  (19)  is  obtained  by  the  image  method  in  conjunction 
with  Fourier  transforms.  An  image  system  consisting  of  an  infinite 
row  of  reflected  images  (shown  in  Fig.  3)  is  introduced  to  satisfy  the 
boundary  condition  on  the  solid  vertical  walls.  The  expression  for  the 
nth  element  of  the  curved  wake  for  such  an  image  system,  based  on 
Eq.  (12),  is 


where 


ln  - U.-O 

(<*v)  - / (Ady/ ds)  2 ds  (20) 

n o k=  — oo  47rJ 

nk 

rnk  = + (y+  2k)2  + (zn-02l* 


z 


— 2b  — • 

<H> 

i i 
1 i 

(H) 

i i 
i i 
J L 

(r 

i 

5 

-D— Y 
1 

Dl 

T 

2h 

ITT 



^-VORTEX  RING 


Figure  3.  Image  system  for  satisfying  boundary  conditions  on 
solid  vertical  wails. 
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In  order  to  satisfy  the  boundary  condition  at  the  ventilated  horizon- 
tal walls  an  additional  potential  (($hW  is  required.  Consequently,  the 
interference  potential  for  the  nth  element  can  be  written: 

<*i)  - (21) 

n n L n nj 


Since  ($v)n  " ^m^n  *-s  known,  it  is  only  necessary  to  determine  (0h^n 
find  the  interference  potential  ((0 i)n) - By  substitution,  the  field  equation 
and  boundary  conditions  for  (0ft)n  are 


v^h  = 0 

(22) 

and 

d$Jdy  = 

—d(j>y/d y = 0 at  y = ±1 

(23) 

d&^/dx.  ± (l/Rldt^/dz  ± AF(?2<^|1/6x5z 

= -[(?</> y '^x  ± (l/'R)«9c6v/(?z  ± \Fd2d>  ^dxdz] 

at  z ± A 
(24) 

<j>h  = 0 at  x i 

- — oe  ; d<j>^/dx.  =*  0 at  x = +« 

(25) 

where  the  subscript  n has  been  suppressed  for  convenience. 

In  order  to  apply  Fourier  transforms  to  Eqs.  (22)  through  (25),  the 
dependent  variable  should  be  absolutely  integrable  (viz,  = 0 at  x = ±®) 

over  the  range  -®  £ x (Ref.  14,  p.  27).  Since  it  is  not,  must  be 

replaced  by  the  axial  perturbation  velocity  caused  by  the  horizontal  walls 
(u^).  From  the  definition, 

d<t>h 
uh  = 17 

(26) 

and  the  inverse  relation 

X — x 

- f uhdx 

Xta-oo 

(27) 

the  field  equation  yields 

V2uh  = 0 

(28) 

with  the  boundary  conditions 

dn^/dy  = 0 at  y * +1 

(29) 

/ X=x  . 1 

uh  ± AFdu^/dz  ± (l/R)d/<JzM  u^dx 

\j(=— » 

^ = -[d(f>v'fdx  ± XFd26v/dxdz.  ± ( 1 /R)(96 y/d z] 

at  z = ±A 
(30) 
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and 


= 0 at  x = too 

Now  it  is  possible  to  use  Fourier  transform  techniques. 

Applying  a complex  Fourier  transform  on  x to  Eq.  (30)  yields 


(31) 


uh  ± AF 


£ * -**>&]  (32> 


where  the  barred  functions  indicate  the  transformed  variables  defined 

by 


g(q.  y,  z) 


g(x,  y,  z)eiqxdx 


(33) 


The  last  term  on  the  left-hand  side  of  Eq.  (32)  reduces  to  (Ref.  15,  pp. 
375  to  376) 

l _J_ 

R (—  iq)  dz 

Now  applying  a finite  Fourier  cosine  transform  on  y and  using  Eq.  (29), 
the  transformed  boundary  condition  at  the  ventilated  horizontal  walls 
can  be  written 


-iquh  ± 


at  z = ±A 


(34) 


where 


g(q,m,z)  = _L_  J g(q,  y,  z)  cos  (nwy)dy 
V 2tt  o 


(35) 


Applying  a complex  Fourier  transform  on  x to  Eq.  (20)  and  using 
the  appropriate  convolution  theorem  yield 


*v(q.  Y.  z)  - j=  / (z.o/  2 — — \ e ' ^qds 

-v  yj*n  i 2n  * \k=-o«  Pk  / 

where  is  the  modified  Bessel  function  of  order  one  and 

= \/(y  + 2k) 2 + (z-O2 


(36) 


(37) 
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Now  applying  a finite  Fourier  cosine  transform  to  Eq.  (36)  gives 


<£v(q,  m,  z)  = / 4>  (q,  y,  z)  cos  (mry)dy  = j-  Ad^/ds  e1[^ds 

~ o y 2^  o " 


r 


where 


oo  1 KjCqpr.) 

H = t 2 / - cos  (m7jy)dy 


k=  — OO  O 


P k 


By  following  the  technique  of  Acum  (Ref.  16),  H can  be  written 


H ' 


where 


Consequently, 


t 4,2  2 2 

I = y q -t-  m 77 


siv(q,  m,  z)  = ~ /-\dy/d3  eH^IVrfds 

V 2w  o 4 | z - £ | 


(38) 


Transforming  the  Laplace  equation  [Eq.  (28)]  yields  the  ordinary 
differential  equation: 


d2u 


d z 


h -2= 

r-  f uk 


(39) 


which  has  the  general  solution: 

uj,  = A cosh  (fz)  + B sinh  (fz)  (40) 

Substituting  Eqs.  (38)  and  (40)  into  Eq.  (33)  gives 

uh  = iq  J"  — Y~  J( A 1 -4-  i A 2)  sinh  (£f)  cosh  (zf)  + (Bj  + B2)  cosh  (^f)  sinh  (zOjeiq^ds  (41) 

where  Aj.,  A2,  Bi,  and  B2  are  given  in  Appendix  B.  Applying  inverse 
Fourier  transforms  on  x and  y yields 

Uh(x,  y,  z)  = -i=  2 j / iqj/  — [(A)  + iA2)  sinh  (CO  cosh  (zf)  (49) 

V 2n  m=o  —00  lo  z ^ ' 

+ (Bj  - iB2)  cosh  (£0  sinh  (zf)]e-it^x“^dsj  cos  (mTry)dq 
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1,  m = 0 

2,  0 


Hence  from  Eq.  (27), 

^ = jk  Jo  j L {o  (^T^)[(A1  + iA2>  8inh  W COsh  (zf> 

+ (Bj  + iB2)  cosh  (£f)  sinh  (zf)Je“i<^x“^ds^  cos  (m^y)dq 
« vfe  J.  j £{/(^)  [(A,  * 1A2)  .ink  W cosh  Ul> 

4 (Bj  4 iB2)  cosh  (£f)  sinh  {zf)j e->q(x-f  cos  (mjry)dq 


(43) 


The  last  term  on  the  right-hand  side  can  be  shown  to  vanish  (Ref.  16, 
p.  61)  as  long  as  the  wake  is  not  horizontal.  For  a horizontal  wake,  it 
is  obvious  that  £ = s,  ? = 0,  and  £ -*■  - thus  the  last  term  of  Eq.  (43) 
becomes 


-lim  J_ 
* " “ y/l  n 


sinh  (z0  cos  (mny)  — 

q 


(44) 


when  the  divergent  improper  integrals  are  evaluated  by  the  Cesaro 
method  (Ref.  17,  p.  361).  Finally,  Eq.  (44)  can  be  written  (Ref.  16, 
p.  61) 

_Z_  £ jBj  sinh  (zf)  cos  (mffv)  AdKds>\  (45) 

4-  0 \'  2it  m=o  1 \ 2 ) 

The  appearance  of  Eq.  (45)  is  a consequence  of  the  absolute  integrability 
requirement  for  the  Fourier  transform.  In  a non-horizontal  wake,  the 
jet  trajectory  is  assumed  to  terminate  at  the  lower  tunnel  boundary  at 
some  finite  downstream  distance,  hence  the  potential  induced  by  the 
horizontal  walls  vanishes  at  infinity  thereby  satisfying  the  absolute  in- 
tegrability requirement.  However,  for  a horizontal  wake,  the  vortex 
cylinder  representing  the  model  trails  to  infinity,  hence  the  potential 
does  not  vanish.  As  discussed  in  Ref.  18,  the  omission  of  the  absolute 
integrability  requirement  on  the  Fourier  transforms  requires  the  addi- 
tion of  a constant  upwash  to  provide  undisturbed  flow  at  x = -“  for  the 
case  of  a horizontal  wake. 
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Finally,  for  a non- horizontal  wake,  the  total  interference  potential 
is  determined  by  combining  the  appropriate  terms  into  Eq.  (21)  yielding 


iA2)  sinh  (£f)  cosh  (znf) 


(Bj  + iB2)  cosh  (£f)  sinh  (fz)j  e-it^x-^ds  j cos  (mffyldq 
‘-^ds 

O \ * / k=-»  2 77  fd 

k#o 


(46) 


For  a horizontal  wake,  the  interference  potential  is  identical  to  Eq.  (53) 
of  Ref.  18.  Unfortunately,  although  a separate  solution  exists  for  an 
exactly  horizontal  wake,  the  general  solution  [Eq.  (46)]  involves  non- 
uniformly  convergent  integrals.  Hence,  in  any  numerical  computation 
the  convergence  problem  becomes  increasingly  severe  as  the  wake  ap- 
proaches the  horizontal  (i.e.,  |30  -*•  ir/2  or  Um/Uj  -►<»).  Therefore,  solu- 
tions determined  from  Eq.  (46)  should  exclude  nearly  horizontal  wakes. 

Lo  (see  Ref.  1),  in  a similar  analysis,  obtained  an  equation  of  the 
same  form  as  Eq.  (46)  for  a continuously  curved  jet  model  in  an  ideal 
slotted  wall  tunnel.  He  then  computed  interference  factors  by  numer- 
ically integrating  the  interference  potential  over  the  length  of  the  jet. 
This  method  proved  to  be  more  laborious  than  the  present  technique  of 
segmenting  the  jet  trajectory,  consequently  his  numerical  computations 
were  restricted  to  determining  the  interference  at  the  model  location 
(x  = 0).  Results  of  the  two  techniques  are  compared  in  Section  3.  0. 


2.4  INTERFERENCE  FACTORS 


The  upwash  and  streamwise  components  of  the  interference  factor 
for  the  model  are  defined  as  (Ref.  5): 


and 


x C ui 
= T~ 


(47) 


(48) 
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where  w[  = 9^i/3z  is  the  boundary- induced  velocity  in  the  vertical  di- 
rection, positive  downward,  and  uk  = 9$j/dx  is  the  boundary- induced 
interference  velocity  in  the  stream  direction.  The  upwash  interference 
factor  Eq.  (47)  is,  by  definition,  four  times  the  conventional  upwash 
interference  factor  for  a wing  (Ref.  5). 

Differentiating  Eq.  (46)  with  respect  to  z and  integrating  over  each 
of  the  N straight  line  segments,  the  total  upwash  interference  factor  is 


oo  oo 


N (Ady/ds)n  j_2X  „ 
w n=o  (Ady/ds)0  J v m=0  ^ 0 


(A,  sinh  (fzn)  i- 

~ 2-  - d2  • K cos  (<Ixn>  ~ Pn  sin 

t cn  n 

A2  sinh  (fz  ) p - Bi  cosh  (fz  ) r 1 

[p„  cos  (qxn)  + Rn  sin  (qxn)J - [Nn  cos  (q*n)  + M„  sin  (qxn)J 


2 ,2 

C J-  d 
n n 


cosh  (f  z ) r -A 

T~a 2 L*M"  C°S  ^ + N"  S!n  (qXn3)dq 

n n / 


4A 


+ T (D"  ~ 3znGn  + 6zn  COS  Fn  « 3 ^n  En> 


(49) 


where  cn,  dn,  Dn,  En,  Fn,  Gn,  Mn,  Nn,  Pn»  and  Rn  are  defined  in 
Appendix  B. 

Similarly,  the  streamwise  interference  factor  is 


S..  = 


N <Ady/ds)n  2A  ^A,  cosh  (fzn)  p , , n . , n 

i.  <15733;  yr  ' / '(  — j2  Lp»  “■  ^ * R» M 


Ag  cosh  (fz„) 


osn  uzj  p -|  B,  sinh  (fz  ) p 

i ~2 L~Rn  C0S  (qXn}  + Pn  sin  (1xn)]  + j LM«  C°S  (qXn>  ~ Nn  sin  (tlxn^ 

n 

B2  sinh  (fz„)i 


„2  H2 

cn  dn 


sinh  ( f z ) p -A 

2 .2  ~ [Nn  cos  <V„)  * V1n  sin 

cn  + n / 

1 2A  00  p t| 

~ — kE  [ZnXnGn  + cos  Pn  sin  Pn  En  “ <zn  sin  Pn  + xn  cos  Pn>Fn]\ 


(50) 
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3.0  NUMERICAL  RESULTS  AND  DISCUSSION 


3.1  JET  TRAJECTORY 

The  interference  factors  for  the  V/STOL  model  can  be  computed  by 
numerical  integration  of  Eqs.  (49)  and  (50)  once  the  jet  trajectory  and 
jet  strength  are  specified.  In  the  present  analysis,  Margason's  (Ref. 

13)  empirical  jet  trajectory  is  used.  In  nondimensionalized  form,  the 
trajectory  is  given  by 


where  Uj  is  the  initial  jet  velocity,  d0  is  the  initial  diameter,  and  |30  is 
the  initial  jet  angle.  The  effect  of  approximating  the  jet  trajectory  by 
line  segments  is  illustrated  in  Fig.  4 where  the  calculated  interference 
factors  at  the  model  location  (x  = 0)  in  a tunnel  with  ideal  slotted  hori- 
zontal walls  (1/R  = 0)  are  shown  for  various  values  of  N.  The  jet 
strength  was  assumed  to  be  constant  along  the  jet.  It  was  determined 
that,  for  the  jet  velocity  ratios  of  interest,  the  curved  jet  path  could  be 
satisfactorily  represented  by  five  straight  line  segments  (N  = 5). 


NUMBER  OF  STRAIGHT  LINE  SEGMENTS,  N NUMBER  OF  STRAIGHT  LINE  SEGMENTS,  N 


a.  Upwash  interference  factor  b.  Streamwise  interference  factor 

Figure  4.  Convergence  ui  aoymented  straight  line  approximation  to  curved  jet  wake. 


20 


AEOC-TR -74-51 


Typical  solutions  for  the  axial  distribution  of  the  interference  fac- 
tors in  an  ideal  slotted  tunnel2  for  a uniform  strength  jet  are  shown  in 
Fig.  5.  The  x location  of  the  jet  intersection  point  with  the  tunnel  lower 
boundary  is  denoted  in  Fig.  5 by  J.  I.  P.  The  agreement  between  the 
present  analysis  and  the  solution  of  Lo  (Refs.  1 or  7)  formed  by  numer- 
ical integration  along  the  curved  wake  is  seen  to  be  excellent. 


a.  Upwash  interference  factor 

Figure  5.  Typical  axial  distribution  of  interference  factors  for 
a constant  strength  jet  in  an  ideal  slotted  tunnel. 

The  terminology  slotted  (or  perforated)  wall  tunnel  refers  to  the 
characteristics  of  the  horizontal  walls.  In  all  cases,  the  vertical  walls 
are  solid. 
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b.  Streamwise  interference  factor 
Figure  5.  Concluded. 


3.2  JET  STRENGTH 

The  calculations  shown  in  Fig.  5 were  done  assuming  a constant 
strength  jet.  As  long  as  the  initial  jet  velocity  is  sufficiently  high  to 
cause  the  jet  to  intersect  the  lower  tunnel  boundary  before  much  decay 
can  occur,  the  average  values  of  the  up  wash  interference  parameter 
calculated  at  the  model  location  are  not  strongly  affected  by  the  constant 
strength  jet  assumption.  However,  an  obvious  paradox  arises  if  a con- 
stant strength  jet  is  assumed.  For  a constant  initial  jet  angle,  the  inter- 
ference factors  for  a weaker  jet  (corresponding  to  a longer  trajectory 
before  intersecting  the  lower  boundary)  are  larger  in  magnitude  than  for 
a stronger  jet  since  the  constant  strength  doublets  are  integrated  over 
a longer  length.  This  paradox  is  illustrated  in  Fig.  6 where  the  axial 
distributions  of  the  interference  factors  are  compared  for  two  different 
values  of  velocity  ratio,  U^/Uj  = 0.212  and  0.311.  The  calculations 
shown  in  Fig.  6 were  made  assuming  a constant  strength  jet  and  an  ideal 
slotted  tunnel.  Although  the  upwash  interference  factor  at  the  model  lo- 
cation (x  = 0)  seems  plausible  for  the  weaker  jet  (U^/Uj  = 0.311)  when 
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a.  Upwash  interference  factor 

Figure  6.  Comparison  of  interference  factors  for  two  different 
constant  strength  jets  in  an  ideal  slotted  tunnel. 

compared  with  the  stronger  jet  (U^/Uj  = 0.212),  the  interference  associ- 
ated with  the  weaker  jet  reaches  a much  higher  maximum  magnitude 
downstream.  This  is  not  consistent  with  any  physical  interpretation  of 
the  real  flow  field.  Furthermore,  as  shown  in  Fig.  6b,  the  streamwise 
interference  factor  has  a higher  magnitude  for  the  weaker  jet,  even  at 
the  model  location.  Hence,  a better  description  of  the  jet  decay  in  the 
tunnel  may  be  as  important  as  the  curvature  of  the  jet  in  the  evaluation 
of  pitching  moment,  tail  force,  and  blockage  corrections.  Unfortunately, 
a generalized  empirical  relation  for  the  jet  decay  is  not  available  at  this 
time.  The  idealized  mathematical  model  used  in  the  present  analysis  is 
a phenomenological  representation  for  a general  V/STOL  model  with  the 
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b.  Stream  wise  interference  factor 
Figure  6.  Concluded. 


circulation  along  the  jet  path  being  related  to  the  entrainment  velocity 
of  a jet  in  cross  flow.  Platten  and  Keffer  (Ref.  19)  indicate  that  the  en- 
trainment velocity  for  a jet  in  cross  flow  is  composed  of  terms  propor- 
tional to  the  velocity  excess  in  the  jet  (producing  entrainment  across  a 
turbulent  shear  layer)  and  to  the  shear  inflow  induced  by  the  vortices 
developed  within  the  jet.  The  rate  of  change  of  this  entrainment  velocity 
is  related  to  the  cosine’ of  the  trajectory  angle.  Hence,  for  an  approxi- 
mation to  the  decay  characteristics  of  the  jet  in  cross  flow  represented 
by  the  present  mathematical  model,  it  may  be  assumed 


(52) 
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which  is  similar  to  the  decay  characteristics  assumed  by  Kirkpatrick 
(Ref.  8).  Actually,  Eq.  (52)  predicts  an  excessively  large  decay  rate 
since  a vortex  motion  is  known  to  persist  in  the  far  downstream  for  a 
jet  in  cross  flow  (Ref.  19).  However,  by  comparing  the  solution  for  the 
constant  strength  jet  with  the  solution  for  the  cosine  decay  model,  an 
upper  and  lower  limit  of  the  interference  factors  for  a more  realistic 
jet  decay  model  can  be  formed. 

Typical  comparisons  between  a uniform  strength  and  cosine  decay 
wake  model  in  an  ideal  slotted  tunnel  are  shown  in  Fig.  7 for  the  upwash 
interference  factor  and  in  Fig.  8 for  the  streamwise  interference  factor. 


a.  UjUj  = 0.099 

Figure  7.  Comparison  of  the  upwash  interference  factor  in  an  ideal  slotted 
tunnel  for  a constant  strength  and  decaying  strength  jet. 
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b.  UJUj  = 0.212 
Figure  7.  Continued. 

For  a high  jet  velocity  as  shown  in  Figs.  7a  and  8a,  the  jet  wake  inter- 
sects the  lower  tunnel  boundary  very  near  the  model,  hence  there  is 
little  difference  between  the  constant  strength  and  decaying  strength  jet 
models  since  the  latter  has  not  had  sufficient  length  to  decay  appreci- 
ably. However,  for  a high  jet  velocity,  it  should  be  recognized  that  the 
solution  may  not  be  applicable  since  the  high  energy  jet  from  a real 
V/STOL  lifting  system  may  cause  flow  reversal  when  it  impinges  on  the 
wind  tunnel  floor  (particularly  for  a closed  tunnel).  Thus,  even  though 
calculations  for  a high  jet  velocity  minimize  the  effect  of  the  jet  wake 
simulation,  the  applicability  of  such  calculations  are  limited  by  the  cri- 
teria for  the  occurrence  of  flow  breakdown  (see  Ref.  6 for  criteria  in 
closed  tunnels). 
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c.  UJUj  = 0.311 
Figure  7.  Concluded. 


As  the  initial  jet  velocity  decreases,  the  disparity  between  the  con- 
stant strength  and  decaying  strength  jet  increases.  Of  course,  if  the 
jet  strength  decreases  enough,  the  V/STOL  model  behaves  more  like  a 
conventional  aircraft,  and  the  details  of  the  jet  modeling  become  irrele- 
vant since  the  interference  factors  can  be  predicted  by  conventional 
theories.  As  the  initial  jet  velocity  approaches  zero,  the  paradox  which 
arises  by  assuming  the  constant  strength  jet  is  emphasized  since  the 
maximum  value  of  the  interference  factors  continually  increases  for  the 
constant  strength  jet  model  instead  of  vanishing. 
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a.  UjUj  = 0.099 


b.  UjUj  = 0.212 

Figure  8.  Comparison  of  the  Streamwise  interference  factor  in  an  ideal  slotted 
tunnel  for  a constant  strength  and  decaying  strength  jet. 
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c.  UjUj  = 0.311 
Figure  8.  Concluded. 


3.3  OPTIMUM  WALL  CONFIGURATIONS 


For  values  of  the  slot  parameter  (P)  at  which  the  interference  fac- 
tors are  minimized,  the  difference  between  the  solutions  for  the  con- 
stant strength  and  decaying  strength  jets  is  also  minimized  (see  Figs. 

7 and  8).  This  is  an  important  observation  since  it  indicates  that  a wall 
configuration  which  gives  minimum  interference  may  minimize  the  inter- 
ference for  a wide  range  of  wake  representations.  The  present  analysis 
can,  therefore,  provide  a useful  tool  for  designing  minimum  interference 
V/STOL  tunnels.  A typical  application  of  this  utility  is  shown  in  Fig.  9 
where  the  slot  parameter  required  to  yield  zero  upwash  interference  at 
the  model  location  is  shown  for  various  values  of  jet  velocity  ratio.  The 
difference  in  the  required  slot  parameter  determined  by  assuming  a uni- 
form strength  or  decaying  strength  model  is  significant.  It  should  be 
noted,  however,  that  the  slot  parameter  required  to  minimize  the  up- 
wash interference  may  not  necessarily  be  the  same  as  that  required  to 
minimize  the  streamwise  interference  or  streamline  curvature  effects. 
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Figure  9.  Slot  parameter  required  for  zero  upwash  interference 
at  the  model  location  in  an  ideal  slotted  tunnel. 


3.4  APPLICATION  OF  THE  QUASI-LINEAR  SLOTTED-WALL 

BOUNDARY  CONDITION 

The  applicability  of  any  theory  is  predicated  on  its  ability  to  simu- 
late the  actual  physical  phenomena.  Binion  (in  Ref.  1)  compared  Lo's 
theory  with  the  experimentally  determined  upwash  interference  on  a 
jet- in-fuselage  V/STOL  model  in  a wind  tunnel  with  slotted  horizontal 
walls  and  solid  vertical  walls.  As  indicated  in  Ref.  1,  the  lack  of  agree- 
ment between  theory  and  experiment  is  associated  with  the  slotted-wall 
boundary  conditions  rather  than  the  theoretical  model.  The  basic  assump- 
tions inherent  in  the  classical  ideal  slotted-wall  boundary  condition  has 
been  examined  in  Appendix  A,  and  a quasi-linear  homogeneous  boundary 
condition  has  been  formulated  to  apply  to  high-lift  V/STOL  models  in 
wind  tunnels  with  ideal  slotted  wall.  The  results  of  the  calculations 
using  the  quasi-linear  boundary  condition  are  compared  in  Fig.  10  with 
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Binion's  experimental  data  and  a theoretical  calculation  using  the  classic 
ideal  slotted-wall  boundary  condition.  In  both  calculations,  the  jet 
strength  was  assumed  constant  since  the  initial  jet  velocity  ratio  was 
sufficiently  high  to  intersect  the  lower  tunnel  boundary  before  much  de- 
cay could  occur.  The  agreement  between  the  theory  with  the  quasi- 
linear  boundary  condition  and  the  experimental  data  is  seen  to  be  ex- 
cellent indicating  that  the  nonlinear  cross  flow  effects  must  be  consid- 
ered in  interference  calculations  in  the  V/STOL  case  and  the  quasi- 
linear  approximation  to  these  effects  derived  in  Appendix  A gives  a 
reasonable  assessment  of  the  high  lift  phenomena. 

THEORY 

U„/Uj  > 0.202  (fe/b  » 0.1015 

h/b  » 0.667  Pq'O* 

BOUNDARY  CONDITION 

CONVENTIONAL,  EQ.  (3),  (l/R-0) 

QUASI -LINEAR,  EQ.  ( 6 ),  (l/R*  - ton  12*) 

EXPERIMENT 

o LU/Uj-  0.202  do/b- 0.1015 

h /b  ■ 0.667 


Figure  10.  Comparison  of  theoretical  solutions  with  experiment. 
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Because  of  the  numerous  tunnel  and  model  parameters  [h/b,  R,  K, 
U./Uj,  dQ/b,  j30,  dr/ds(s)],  it  is  not  feasible  to  provide  a comprehen- 
sive set  of  charts  or  tables  of  the  interference  factors.  However,  the 
calculations  can  be  performed  quickly  for  any  given  set  of  variables  and 
a listing  of  the  Fortran  IV  computer  program  for  calculating  the  inter- 
ference factors  is  included  in  Appendix  C to  facilitate  usage  by  inter- 
ested readers. 


4.0  CONCLUSIONS 


A theoretical  investigation  of  wind  tunnel  wall  interference  on 
V/STOL  models  in  rectangular  ventilated  wind  tunnels  has  resulted  in 
the  following  conclusions: 

1.  The  assumption  of  a constant  jet  strength  produces  a 
paradox  in  that  the  maximum  value  of  the  interference 
factors  increases  as  the  initial  jet  velocity  decreases. 

2.  A detailed  description  of  the  V/STOL  model  jet  decay 
may  be  as  important  as  the  curvature  of  the  jet  in  the 
evaluation  of  pitching  moment,  tail  force,  and  blockage 
corrections.  Furthermore,  some  approximation  to  the 
jet  decay  is  required  to  overcome  the  paradox  which 
arises  by  assuming  the  constant  strength  jet. 

3.  A wall  porosity  configuration  which  gives  minimum 
calculated  interference  may  minimize  the  interference 
for  a wide  range  of  wake  representations.  Conse- 
quently, the  details  of  the  theoretical  V/STOL  model 
are  not  essential  to  determine  the  porosity  criteria 
for  a minimum  interference  tunnel. 

4.  The  nonlinear  cross-flow  velocity  heretofore  not 
treated  in  wind  tunnel  interference  calculations  is 
of  importance  in  the  analysis  of  high-lift  V/STOL 
models.  A quasi- linear  homogeneous  boundary 
condition  for  ideal  slotted  walls  has  been  formu- 
lated to  account  for  the  cross-flow  effects  in 
V/STOL  testing,  and  application  of  the  new  boundary 
condition  has  yielded  a significant  inprovement  in 
the  correlation  between  theory  and  experiment. 
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APPENDIX  A 

DEVELOPMENT  OF  QUASI-LINEAR  SLOTTED-WALL  BOUNDARY  CONDITION 

An  equivalent  homogeneous  boundary  condition  for  slotted  walls 
[Eq.  (5)]  has  been  derived  by  various  authors  (e.g.,  Refs.  10  and  20). 

In  each  case,  it  was  assumed  that  both  the  slot  flow  velocities  and  the 
perturbations  from  the  mean  flow  are  small  compared  with  the  undis- 
turbed tunnel  velocity.  However,  for  high-lift  V/STOL  models,  cross- 
flow  velocities  in  the  test  section  may  become  appreciable  with  the  ve- 
locity near  the  slots  being  quite  high.  Therefore,  near  the  slotted  walls, 
the  quadratic  terms  of  the  Bernoulli  equation  may  no  longer  be  neglected 
as  in  the  derivation  of  Eq.  (5).  The  derivation  given  below  heuristically 
makes  an  approximation  to  the  quadratic  terms  of  the  Bernoulli  equation. 

A fundamental  assumption  is  made  that  the  higher  cross-flow  ve- 
locities near  a slotted  wall  induced  by  a lift- augmented  model  may  be 
approximated  by  a uniform,  constant  cross-flow  related  to  the  increased 
equivalent  circulation  induced  by  the  lift  augmentation  device.  High- 
lift  V/STOL  models  characteristically  have  a large  zero-lift  angle  (a0) 
caused  by  any  lift  augmentation  device.  It  is  assumed  that  a high-lift 
model  can  be  represented  as  a conventional  wing  at  zero  angle  of  attack 
embedded  in  a stream  with  a uniform  cross-flow  whose  velocity  is  re- 
lated to  the  zero  lift  angle  as  shown  in  Fig.  A-l. 

HIGH  LIFT  WING  IN  FREE  - AIR 


CHORD 


CONVENTIONAL  WING  IN  CROSS  FLOW 
APPROXIMATION  TO  HIGH  LIFT  WING 


V8 


V,=  U„  lin  do 


Figure  A-1 . Representation  of  a high-lift  wing  as  a conventional 
wing  in  uniform  cross  flow. 
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When  this  equivalent  system  is  placed  in  the  wind  tunnel,  the  free- 
stream  static  pressure  far  upstream  is  given  as 


P=o  = P- 


£ (U2  cos2a  + V2) 


= P. 


- £ u2 


(A-l) 


The  static  pressure  at  a point  inside  the  test  section  at  the  wall  is  given 

by 

P = Pt^  “ | [(U^  c°s  a0  + u)2  + (Ve  + w)2  - v2] 

(A-2) 

= Pi  - ^ [u*  + cos  a0  + 2wVe  + u2  + v2  + w2] 

Now,  if  it  is  assumed  that  the  pressure  just  outside  the  tunnel  wall  is 
maintained  at  p^,  the  pressure  difference  across  the  wall  is  obtained 
from  Eqs.  (A-l)  and  (A-2) 

Ap  = p - = ~2.[2uUeo  cosaD  + 2wVe  + u2- v2  - w2]  (A- 3) 


Further,  if  it  is  assumed  that  u,  v,  w are  much  smaller  than  U^,  then 

Ap  = cos  aQ  + wVe)  (A- 4) 

Obviously,  for  a conventional  lifting  model,  Eq.  (A- 4)  reduces  to  the 
well  known  result  for  small  disturbance  theory. 

In  order  to  determine  the  required  boundary  condition  at  a slotted 
wall,  it  is  necessary  to  find  an  expression  which  relates  the  pressure 
difference  across  the  wall  to  the  flow  through  the  wall.  Consider  a thin 
slotted  wall  in  a field  of  flow  with  a uniform  velocity  normal  to  the  wall 
in  a transverse  plane  at  -ZQ  (see  Fig.  A-2).  Because  the  flow  pattern 
is  the  same  for  each  slot  it  is  permissible  to  study  a single  channel  such 
as  the  one  in  which  approximate  streamlines  have  been  sketched  in  Fig. 
A-2.  The  analysis  of  the  flow  near  the  slotted  wall  presented  here  fol- 
lows closely  the  development  of  Davis  and  Moore  (Ref.  20).  The  kinetic 
energy  enclosed  in  a region  of  the  flow  bounded  by  the  transverse  plane 
at  -ZQ,  by  the  "walls"  of  the  channel,  and  by  the  slot  is 


Kinetic  Energy 


Pff, 


A AX 

<t>  ■=—  dA 

On 


(A-5) 
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Figure  A-2.  Schematic  of  flow  field  perpendicular  to  slotted  wall. 


where  the  area  of  integration  (A)  consists  of  a surface  of  unit  depth  nor- 
mal to  the  plane  of  the  page  encircling  the  region  shown  in  Fig.  A-2. 
Since  dQf  9n  = 0 at  the  channel  walls  and  ♦ = 0 at  the  slot,  these  regions 
contribute  nothing  to  the  integral.  Furthermore,  if  the  transverse 
plane  at  -ZQ  is  sufficiently  far  away  from  the  slot,  the  potential  $(-Z0) 
is  essentially  constant  in  this  plane,  thus 


Kinetic  Energy  = H p d>(-Z0)V'nA 


(A-6) 


where 


v ■ a 


dA 


is  the  volume  flow  rate.  The  potential  at  -Z0  in  the  presence  of  a 
slotted,  wall  has  been  derived  by  Chen  and  Mears  (Ref.  11)  by  replacing 
the  slotted  wall  by  an  infinite  series  of  doublets.  The  result  is 


<D(_Z0)  = Vn(-ZC*K)  (A- 7) 

where  K is  the  geometric  slot  parameter  given  in  Eq.  (A-4).  Therefore, 
Eq.  (A-6)  can  be  written  as 


Kinetic  Energy  = ^p(-Z0  + K)V^A 


(A-8) 
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The  portion  of  the  total  kinetic  energy  regarded  as  being  caused  by  the 
presence  of  the  slotted  wall  is  pKV^/2  per  unit  area. 

Now  if  the  slotted  wall  is  replaced  by  an  equivalent  homogeneous 
wall.  in. the  same  manner  as  originally  proposed  by  Busemann  (Ref.  20), 
the  energy  per  unit  area  at  the  wall  is  pKVn/2.  The  pressure  differ- 
ence across  the  slotted  wall,  which  acts  in  a direction  normal  to  the 
wall  surface  must,  in  potential  flow,  be  equal  to  the  rate  of  change  of 
momentum  associated  with  the  presence  of  the  slots,  or 

n d v_ 

- i,  W - 1*  sr  (A'9) 


If  it  is  recognized  that  Vn  = Ve  + w,  the  derivative  is  given  as 


Dt 


Dw  dw  /rI 

157=  a7+  (u» 


COSao  + u)| 


(A- 10) 


For  steady  flow,  and  to  the  same  order  of  approximation  used  in  deriv- 
ing Eq.  (A-4),  the  pressure  difference  across  the  wall  is  related  to  the 
velocity  through  the  wall  by 

Ap  = pK^  cosa0|j  + Ve  (A- 11) 


Equating  Eqs.  (A- 11)  and  (A-4)  yields 


u + 


Ve 


U cos  a 


w + 


,,  dv,- 

Kax 


KVe 
U cos 

oo  O 


(A- 12) 


or  in  terms  of  the  perturbation  potential: 

do  Ve  d<I>  (920  KVc  d20 

dX  + U cos  a dZ  +KdXdZ  + 0 cos  a 

oo  o oc  o O L 


(A- 13) 


Equation  (A- 13)  is  a linear  homogeneous  boundary  condition  for  the 
slotted  wall  in  the  presence  of  a uniform  cross  flow  with  velocity  Ve. 
By  applying  the  chain  rule  along  a streamline,  it  is  seen  that 

92<t>  _ <5w  dw  dX 
dz2  dZ  = dX  dZ 


However,  by  irrotationality,  9w/9X  = 9u/9Z.  Thus, 

920  _ 9X  _ 3ZQ  dX_ 
dz2  ~ dZ  dZ  ~ dZdX  dZ 
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Therefore,  Eq.  (A- 13)  becomes 

<90  Ve  ao  Ve  dx\  a2o 

ax  + u cos  a az  + \ + u coso.  az/axan 

00  o oo  o 


(A- 14) 


Furthermore,  if  it  is  assumed  that  the  streamlines  at  the  slotted  wall 
are  controlled  principally  by  the  cross-flow  velocity  (Ve),  then  the  angle 
of  the  streamlines  along  the  X- direction  is  approximately  constant  at 
the  wall  and  equal  to  aQ.  Therefore, 


U 


oo 


Ve 

cos  aD 


ax  _ u~  sincto 

az  U cos  a_ 

oo  O 


cot  aD 


1 


(A- 15) 


Thus,  the  slotted-wall  boundary  condition  for  a wall  placed  normal  to 
the  direction  of  lift  takes  the  form  of  the  classic  slotted/ porous  boundary 
condition  derived  by  Baldwin,  et  al.  (Ref.  10): 


ao  j_  ao  „ a2o  _ 
ax  + Re  az  + e ,jxaz 


(A- 16) 


where  the  pseudo  porosity  parameter  (1/Re)  is 


J Ve 

Re  = UM  cosa0 


tan  a0 


and  the  effective  slot  parameter  is 


(A-17) 


Ke  = 2K 


(A- 18) 


Equation  (A- 16)  is  the  slotted-wall  boundary  condition  which  can  be  used 
to  determine  the  interference  effects  on  V/STOL  models  with  high  cross- 
flow  velocities. 
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APPENDIX  B 
LIST  OF  FUNCTIONS 


The  expressions  for  the  functions  Ai,  A2,  Bi,  and  B2  in  Eq.  (41) 


are: 


Ai  = 


[(l  - AfF)(coah  (Af)  + AFf  sinh(Af))-  (f/gR)2  sinh  (Ao]e'Xf 
[(cosh  (Af)  + AFf  sinh  (Af))2  + (f/qR)2  sinh2  (Af)] 


A2  - 


-(f/qR) 


[(cosh  (AO  - AFf  sinh(Af))2+  (f/qR)2  sinh2  (Af)J 


(B-l) 

(B-2) 


B = [(1  - AFfXsinh  (Af)  + AFf  cosh  (Afj)-  (f/qR)2  cosh  (Af)]e'Xf 

1 [(sinh  (Af)  + AFf  cosh  (An)2  + (f/qR)2  cosh2  (Af) 


B, 


(f/qR) 


(sinh  (Af)  + AFf  co  sh  (Af))2  + (f/qR)2  cosh2  (Af)] 


(B-3) 

(B-4) 


The  expressions  cn,  dn,  Dn,  En,  Fn,  Gn,  Mn,  Nn,  Pn,  and  Rn  in 
Eqs.  (49)  and  (50)  are 

Cn  = f cos  /3n 
dn  = qsir>0n 


Dn  = (4a 


2 ^n  + bn  _ K_ 

an  " bn)  _ '"o 


(B-5) 

(B-6) 

(B-7) 


] 

bn 

(4an  + 

bnZ) 

<2En+V| 

, 3 
no 

2 12 

(4an- 

bn> 

12  J 

— | 

bn 

1 + 

4wn  + bn 

2 

(4an~bn> 

■ 12,”* 

4 n - b2 
n n 

“ 3 

«4i‘n-bn2) 

Van 


F 


" 3r  3 
n _ 


1 + 


bn<2t'n  + bn> 


(4a  - b‘) 

n n 


8b  (2£  + b) 
n n n 


3(4an"bn)  rn„ 


3a 


3/2 


(4a.  - b‘) 


(B-8) 


8 b, 


2 


nj  3(4an-bn)Van 

(B-9) 
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( 

it 


(2f  r i>J 


! (l.i  - li‘-)r  ’1 
n n n 


i 


<■».«  - 1.*) 
n n 


8-'n 


(4  it 


n 


(B-10) 


M. 


cn  Sinl1  (cn  (’n)  sin  (d„‘V  “ dtt  tnsh  <cnP„>  tos 


(B-ll) 


\l  " ‘it  cos  sinl'  ! d„  c,,sI’  sin  <dnrn) 

P„  - d„  sinl.  (•:/„)  L-..H(dnf„)  - cn  c..-*  (cnFn)  fin  (anfn) 

Hn  ■ "S.  1 ‘’n  si,,h  (t./J  *iin  (d,.(V  * cn  Cosl>  <c„fn>  LOS  (<1nFn) 

where 

2 .1  2 2 
an  " x..  * lk  + zt. 

[,n  = -2{vn  sin  Pn  ~ 'n  C0R  PJ 

and 


V«„  - bBr„ 

i n n n 


(B-12) 
(B-13) 
(B- 14) 

(B- 15) 
(B- 16) 

(B-17) 
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APPENDIX  C 

FORTRAN  IV  PROGRAM  FOR  THE  IBM  370/155 


I.  DESCRIPTION  OF  PROGRAM 

A.  MAIN  PROGRAM 

The  main  program  receives  the  input  data,  controls  the  calculation 
of  the  jet  trajectory  and  the  interference  factors,  and  prints  the  trajec- 
tory location  and  interference  factors  as  output  data. 

B.  EVAL 

This  subroutine  calculates  the  coefficients  Dn,  En,  Fn,  and  Gn  and 
forms  the  terms  of  the  k -series  of  Eq.  (49)  or  Eq.  (50). 

C.  FCT 

This  FORTRAN  Function  calculates  Ai,  A2,  Bi,  B2,  cn,  dn,  Mn, 

Nn,  pn>  and  Rn  and  forms  the  integrand  for  the  inverse  Fourier  integrals 
of  Eq.  (49)  or  Eq.  (50). 

D.  TF 

This  FORTRAN  Function  calculates  the  trajectory.  For  the  pres- 
ent analysis,  Margason's  trajectory  [Eq.  (51)7  was  used.  Any  trajec- 
tory equation  of  the  form  X = X(Z,  U^/Uj,  dQ/b,  j3Q)  can  be  substituted 
into  this  function  by  the  User. 

E.  BETAF 

This  FORTRAN  Function  calculates  the  local  slope  of  the  wake 
trajectory.  The  equation  used  must  be  compatible  with  the  Function  TF. 

F.  DF 

This  FORTRAN  Function  calculates  the  decay  of  the  jet  along  the 
jet  wake.  Any  user  specified  function  of  the  form  DF  = DF(TF,  BETAF) 
can  be  used. 
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G.  LAG  32 

This  subroutine  performs  the  numerical  integration  of  the  inverse 
Fourier  transforms  of  Eq.  (49)  or  Eq.  (50)  by  a standard  32 -point 
Laguerre  integration  scheme. 

II.  INSTRUCTIONS  FOR  USE 

A.  DEFINITION  OF  INPUT  VARIABLES 


N 

Number  of  straight  line  segments 
mating  the  curved  trajectory. 

approxi- 

XLAM 

Tunnel  height-to-width  ratio,  X = 

h/b. 

UE 

Jet  velocity  ratio,  U^/Uj 

DO 

Ratio  of  initial  jet  diameter  to  tunnel  semi- 
height, dQ/b 

BO 

Initial  jet  angle,  J3Q>  deg 

NM 

Number  of  x locations 

XMIN 

Initial  x/b  location 

XDEL 

x/b  increment 

P 

Slot  parameter 

Q 

Porosity  parameter 

B.  ORDER  OF  DATA  DECK 


Columns 


Variable 


Card  1 

0-5 

N 

6-15 

XLAM 

16-25 

UE 

26-35 

DO 

36-45 

BO 

Card  2 

0-5 

NM 

6-15 

XMIN 

16-25 

XDEL 
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The  input  values  for  P and  Q are  internally  generated  in  the  main 
program.  They  are  assigned  the  values  0.001,  0.2,  0.4,  0.6,  0.8,  and 
1.0.  Particular  values  for  P and  Q can  be  generated  by  appropriately 
modifying  the  Data  AH/ card  (line  7 of  the  main  program). 

C.  DATA  VALUES 

The  variables  N and  NM  must  be  right  justified  in  their  fields,  and 
punched  without  a decimal  point.  The  variables  XLAM,  UE,  DO,  BO, 
XMIN,  and  XDEL  must  be  punched  with  the  decimal  point,  but  need  not 
be  right  justified. 

D.  INPUT  AND  OUTPUT  SAMPLES 

A sample  input  deck  is  shown  following  the  program  listing.  The 
output  for  this  case  is  also  shown.  A description  of  the  output  follows: 


Line 

1: 

Input  variables  N,  h/b,  U^/IL,  dQ/b, 
and  /30 

Line 

2: 

Trajectory  description 

Line 

3: 

Jet  decay  description 

Line 

4-10: 

Jet  trajectory  geometry,  j3n,  ?n,  and  l?n  for 
n = 0,  1,  . . . , N 

Line 

11: 

Slot  and  porosity  parameters,  P and  Q, 
respectively 

Line 

12: 

Values  of  6w^»  *^w2'  ^ui»  ^u2*  ^w»  ^u  at 
the  x/b  locations  where  <5wi  and  6Ui  are 

the  contributions  from  the  horizontal  walls, 
6W 2 and  6u2  are  the  contributions  from  the 
vertical  walls,  and  6W  and  6U  are  the  total 
interference  factors. 
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V/STOL  MODEL  WALL  INTERFERENCE 
USING 

SLOTTED/POROUS  BOUNDARY  CONDITICN 


IMPLICIT  RE  AL  *8 1 A— H , 0-  Z ) 

DIMENSION  ZETAI 101 ,XI ( 10) , A( 10)  ,8{ 10 ) , C 1 10 ) , DW1 1 10 ) ,0W2( 101, 

1DUK 10) ,DU2( 10) ,B£TA(10)«SM8(10),DFFI10) 

COMMON  /PASS/  BO. DO.UE.XIN 

COMMON  XI 10),Z(10)  .XL 1 10) . SBCIO ) ,CB 1 10) ,F,R IN.XLAM. M. I. PI, PI  2 
DIMENSION  AH ( 6 ) 

DATA  AH/O. 0010+0, 0.20+0, 0.40+0, 0.60+0, 0.80+0, 1.00+0/ 

CALL  ERRSFT(207,256,-l,l) 

CALL  ERRSF  T ( 266  ,2  56, **1,1) 

CALL  ERRSET(208»256,-l,l) 

900  REAO  (5, 500, END-999)  N.XLAM.UE ,DO,BO 
WRITE  (6,130)  N.XLAM.UE, 00, BO 
READ  (5,200)  NM.XMI N.XOEL 
WN=N 

P 1=3. 141 592626535 8979 
P I2=P 1/2. OD+O 
XIN=TF(— XLAM.O) 

CALL  OFN 
ZETAN*-XLAM 
NP1=N+1 
WRITE  16,121) 

DO  12  1=1, NP1 

ZETAI I )=ZETAN*(I-1)/WN 

XI(  I)=TF(ZETA(I),0) 

DFF(  I)  = DF(ZETA(I) , XI  ( I ) ) 

BFTAI  I)  = BFTAF(ZETA(D) 

BETA  ID*  80/57. 29577951310+0 
SB(  1 ) =DSI N( BETAf I ) ) 

CB(  I ) =DCOS( BETA( I ) ) 

12  WRITE  (6,120)  BETA(I) , XI ( I ) , ZETAI II 
TEMP- ZE  TAN/WN 

DO  13  I =2 ,NP1 

T=X  I ( [l-XIIl-U 

SMB ( I )=OSQRT(TEMP*TEMP+T*T) 

T=DATAN(-TEMP/T) 

B(I  ) = PI2-BETA(I-1)-T 
C( I ) =-PI2+BETA( I )+T 

13  A(I)=PI-BETA(I)+8ETA(I-1) 

XL< 1)=SMB(2)*0SINIC(2) )/DSlN(A(2) ) 

DO  14  1=2, N 

14  XL(I)  = SMBm*DSIN(B(I))/DSIN(A(I)!  + SMB(I  + l)*DSIN(C(  1+1 ) l/DS  INI  All  ♦ 
ID) 

XL(NPI)=SMR(NP1)*0SIN(B(NP1) )/DSI NC A (NPl) ) 

DO  998  KR=1 ,6 
P=AH(KR) 

DO  998  KS=l,6 
(J=AH(  KS ) 

WRITE  (6,1561  P,0 
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F=l.  00+0/ P-1. OD+O 
KIN=1.0D+0/Q-1.00+0 

XL(  NP1)  =SMB  INP1 1 *05 1 N(  HINPlt  I/0SIN(A(NP1I) 

WRITF  (6,112) 
no  60  J=1 , NM 
XX=XMIN  + ( J-l ) *XDEL 
X(l)=XX 

z(  n=o.oo+o 

00  16  1=2, NP1 

IF  ( I .NE.2)  GO  TO  15 

X(I J=Xl( I - 1) *SB( I -1 ) 

Z(  i j=xl(  i-h*cb(  i — i » 

GO  TO  16 

15  xm  = x<  I-lJ+XLt  l-l)*SB(I-l) 

Z(II=Z( I-1)+XL(  1-1)  *C8«  1-11 

16  CONTINUE 

00  17  1=2, NP1 

17  X(!)=XX-X(I) 

1=1 

is  nuim=o. oo+o 

M=0 

19  CALL  LAG32C  0.0D+0,Y,1) 

190  IF  (M.NE.O)  Y=2.00+0*Y 

OUll  I)=0>J1(1)+Y 
IF  (M.GT.250)  GO  TO  20 
IF  CM.LT.3)  GO  TO  80 

IF  (OABS(Y) • LE . 1. 00— 09+DABS ( OU 1(1)1)  GO  TO  20 

80  M=M  + 1 

GO  TO  19 

20  ouim=2.on+o*xLAM*ouim/Pi 

IF  ( I .E0.NP1)  GO  TO  83 

1 = 1 + 1 

GO  TO  18 
33  1*1 

21  DU1( 11=0.00+0 
M=0 

22  CALL  LAG32(0.00+0,Y,2) 

IF  IM.NE.OJ  Y=2. OO+O+Y 
DUl ( I )=0W1( I )♦ Y 

IF  CM.GT.250)  GO  TO  23 
IF  (M.LT.3J  GO  TO  81 

IF  (DABS(Y).LE.l. 00-05*0 ABS < DW 1(1))}  GO  TO  23 

81  M=M+ 1 

GO  TO  22 

23  0M1(  I )=2.0D+0*XLAM*DH1(  IJ/PI 
IF  ( I.FQ.NPl)  GO  TO  2* 

1 = 1 + 1 
GO  TO  21 

24  1 = 1 

93  DW2 ( I)=0.0 
K=1 

25  CALL  E VAL  ( K , 1 , AIMS  ) 

0W2 ( I)=0W2(I)+ANS 

IF  (K.LE.3)  GO  TO  51 
IF  (K.GT.250)  GO  TO  26 

IF  (DABS(ANS).LF. l.0D-05*DABS(DW2(I)> ) GO  TO  26 
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51  K*K+i 
GO  TO  25 

26  DW2(I  )=4.00+0*XLAM*DW2(I)/PI 
IF  ( I .EQ.NPl)  GO  TO  27 
1*1*1 

GO  TO  93 

27  1*1 

94  DU2(  I )=0.0D+0 
K=1 

28  CALL  FVAL1K.2.ANS) 

L)U2(  I I =DU2(  I l+ANS 

IF  (K..LE.31  GO  TO  50 
IF  (K  .GT.250)  GO  TO  29 

IF  (DABS! ANSI.LE. 1. 00-05*DABS( DU2  ( I ) ) ) GO  TO  29 
50  K=K+ 1 
GO  TO  28 

29  OU2(II*-l2.0*DU2t!l/PI 
IF  ( I .EQ.NP1I  GO  TO  30 
1*1  + 1 

GO  TO  94 

30  SUM1=0. 00+0 
SUM2*0. 00+0 
SUM  3=0. OD+O 
SU**4-0.0D+0 
on  40  1*1. NP1 

SUM  1=  SUM 1 + OFF ( 1 1+0W1 ( I ) 

SUM2* SUM2+DFF ( I » *DM2C 1 1 
SUM 3= SUM3+0FF ( I ) *DU1 ( I I 
40  SUM4*  SUM4+DFF  Cl) *DU2 (II 
0W=SUM1+SUM2 
OU* SUM3+SUM4 

WRITF  (6.120)  XX.SUMl,SUM2fSUM3.SUM4f0W.0U 
60  CONTINUE 

998  CONTINUE 
GO  TO  900 

999  STOP 

101  F0RMAT(///,5X, 'TRAJECTORY  DOES  NOT  INTERSECT  WALL'.//) 

102  FORMAT! ///,5X, 'TRAJECTORY  INTERSECTS  WALL'.//) 

156  FORMAT(/,5X.»P=',F5.3,5X,'Q*',F5.3,/) 

112  FORMA T(//,17X,'X* ,15X , 'OWl • , 15X , 'DW2 ' , 1 5X , ♦ 0U1' , 15X, • 0U2' , 16X, 

1 • DW ' , 16X, 'OU',/) 

122  FORMAT! //  ,24X.'A'  ,17X , *8 • , 17X , • C' • 1 5X, • SMB' ) 

121  FORMAT! 21X, 'BETA'.lbX.'Xl*. 14X  » 'ZET  A* ) 

120  FORMAT! 7X.7D18. 7) 

130  FORMAT  ( 1H 1 . ' N*» ,1 5 , 3X, ■ LAMBDA*' ,F 10.5 . 3X , 'UE=* . F10.5, 3X, » DO-' , 
1 F10.5,3X,'B0*’,F10.5) 

200  FORMAT  (T5.2F10.0) 

' 500  FORMAT  ( I 5.4F10.0) 

END 
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subroutine  EVALU.ICASE,  ANS) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON  XN(IO) ,ZN(lO) ,XL(10) ,SB ( 101 . CBC 10 ) ,F,R IN.XLAM, M, I, PI, PI  2 
AN=XN(  I ) *XN(  I >4-4.0040*K*K4-ZN(I  )*ZN(  I) 

SA  = 0 SORT! AN ) 

A32=SA**3. 

BN=-2.0D+0*IXNI I! *SB( IJ-ZNC I )*CB(II  ) 

6S=8N*RN 
V=4.0D+0*AN 
T =V-R  S 
U=V+BS 
8SQA=8N/SA 

RN=  AN4-XL I I I *( BN+XLI I I I 
SR=D  SORT!  RN ) 

R3=SR*SR*SR 
S=2.00+0*XL( l )+BN 
ON=2.004-0*(  S/SR-BSQAJ/T 

FN=  ( ( 1. 004-0+8  S/T I /A32+8.00+  0*BS/(T*T*SA)-(1.0D+0+8N*S/T  )/R3  - 

1 8.00+0*BN*S/(T*T*SRI  1/3.0 

FN*-FN 

GN=2. 00*0*1 S*( 1.00+ 0+8. 0*RN/T)/ R3-BN*( 1 .00+0+8. 0D+0*AN/T )/A321  / 

1 (3.00i-0*T) 

EN*  (-XL ( I )+(BN+U*S/T)/6. 00*0) / ( 2.0D+0+R3) +2 • OD+O*  (U*S/SR-U*  BSUAI/ 
l ( 3 • OD+O+T*  T)-BN*  ( 1. 00*04-11/  T) / ( 12 .9 *A32  ) 

CC=CP(I  ) 

IF  ( ICASE.E0.2)  GO  TO  20 

ANS=0N-3.00+0*ZNl  I ) * ( ZN ( I )*GN+6.0D+0*CC*FN)-3.0* EN*CC*CC 
RETURN 
20  SSsSB(  I ) 

ANS*ZN(  M*XNm*GN-(ZN(!  )*SS4-XN(  I ) *CC  )*FN+CC*SS*EN 

RETURN 

END 
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rUNC  T ION  FCr(X,KI 
[MPLICTT  RFAL*ril A-H.O-Z) 

COMMON  Yl 10» ,Z(10) , XL (10 1 , SB  1 10 1 ,CB(10J ,F ,k IN, XL AM, M, I,PI,P!2 
PM=DS0R  T ( X*X*M*M*PI*Pl  I 
PL=PM*XLAM 
£ PL  =0EXP{ -PL ) 

SPL=QS1  NH(PL) 

CpL=l)CUSH  ( PLl 
PFL=F*PL 
PSI=1. J0*0-PFI 
P1=CPL*PFl*SPL 
P2=  SpL*PF  L*C  dL 
0N=X#  SB ( I I 
CN=PM*f9«  I I 
3X=X*Y<  I I 
SQX=0SI N( OX) 
r.ox=ocos<  ox) 

DL=DN*XL<  I ) 

CL=CN*XL<  I ) 

SPL2=SPL*SPL 

CPL2=CPL*CPL 

SJL=DSIN(DL> 

COL=OCOS( 01 ) 

SHU  =nSlMH(OL) 

CHOL=OCOSH(Ol) 

SCL=0  SI N( CL ) 

CCL=OCnS(CL) 

SHCL=OSINH(CL I 
CHCI.=ncnSH(CL) 

CMN=ON*CN*SHC  l *sdl— dn*chc  l*col 
ENN«CN*CDL*SHCl+DN*CHCL*SDL 
PN=ON*SHCL*CDl-CN*CHCl*S0l 
H.N=-CN4-0N*SHCL*S0L«-CN*CHCL*CDL 
ZPH=PM*Z« I ) 

CHZPM=QCOSH( ZPM) 

SHZPM=OSI NH(ZPM) 

CNDN=CN*CN*DN*DN 

S=PM*RTN/X 

S2=S*S 

UENOM=  Pl*Pl>S2*SPL2 

A 1 1=  (PSI*P1-S  + SPL)*EPL  / OENOM 

A L2=  -S  / DENOM 

f>ENOM=P2*P2+S2*CPL2 

Bll  = <-PSI*p2«-S2*C!>l)*EPl  / DENOM 

312=  S/  DENOM 

IF  (K.E0.2)  GO  TO  20 

T L=A1 1*CHZPM*( PN*COX^RM*SOX) 

T2=A12*CHZPM*( PN*SQX-RN*COX) 

T3=B11*SHZPM*JEMN*CQX-ENN*S0X) 

T4=B12*SH7PM*CENN*CQX«-F.MN*SQX) 

FCT=  (TltT2+T3+T4l *X/CNDN 
RETURN 

20  Tl=All*SHZPM*(PN*SOX-RN*CQX) 

T2=-A 12*SHZpM«( PN*CQX*RN*SOX ) 

T3=Bl  L*CHZPM*(FNN*COX-*-EMN*SQX) 

TA=HL2*CHZPM*<  ENN*SQX-EMN*CQX) 

FCT=  (TH-T24-r3*TA)*PM/CN0N 

RETURN 

FNO 
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FUNCTION  TF<X,K) 

IMPLICIT  RCAL*8(A-H,0-Z) 

COMMON  /PASS/  H0,D0,UF.,XIN 
OATA  11/0/ 

IF  I I l.EQ.O)  WRITS  (6,1001 
I 1*1 

100  FORMAT  ( /// ,5X,  • MARGASOM'S  TRAJECTORY* I 

TF=-0.250*0*UC  *UE  *X**3/ ( OO*DO*DCUS0  (001**2) 
1) 

TF  JK.EQ.O)  RETURN 

TF=TF-XIN 

RETURN 

END 


FUNCTION  RFTAF(X) 

IMPlfCIT  RFAL*8(  A-H,f)-Z) 

COMMON  /PASS/  rt0,00,Ufc 

4 = 0.  7l>0*0*Ur  *UF  *X*X/  ( 00*00  *DCOSDI  00  ) **2  ) 
OF  TAP  = OATANI 4) 

RETURN 
F NO 


FMNC  THIN  DP (X,Y) 

IMPLICIT  RCAL*fi( A-h,0-Z) 

COMMON  /n  A S S/  AU,O0,UF. 

OF* 1 .0  )P0 
RETURN 
ENTRY  JPN 
WRITE  (6,100) 

100  FQRMAK ///, 5X, ‘CONSTANT  STRENGTH  JET‘1 
RFTUPN 
f NO 


- X*OS IND( BO ) /DCUSO( BO 


OS  I NO ( 00 ) /OCOSD( 00 ) 
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SUHRnUTINE  LAG32C ft » V »K> 

DOUBLE  PRECISION  X,Y,A 
FCTLCX)=DEXPIX) *FCT(X*A,K» 

X=. 111751398097938D  03 

Y=.451053619389897D-47*FCT1 CX) 

X=.  98829542  86 82 84 OD  02 

Y=Y+. 13 3861 69421 06260-41 4FCT1CX) 

X=. 8873534041789240  02 

Y=Y4. 267151 121924014D-37*FCTl(X) 

X=. 80187446977913  50  02 

Y=Yt. 1192248760098220-33 *FCT1CX» 

X= • 7268 7o2 8090662  7D  02 

Y=Yf. 19 1337549445422 D-30*FCT1(X» 

X=. 6597537728793510  02 

Y=Y+. 14 18560545463040— 27 *FCT1 (X ) 

X=. 5989250916213400  02 

Y=Y+. 5661 29 413 03 973 60-25 *FCT1JXI 

X=. 5433372133339690  02 

Y=Y+. 1346 982586b3 740 D-22*FCT1CX1 

X=. 492 243 94 9 87 308 6D  02 

Y-Y+. 20544296 7378805 D-20*FCT1(X» 

X=. 4450920793575490  02 

Y=Y  + . 21 19792290163620-1 8 *FCT1 C X ) 

X=. 4014571977153940  02 

Y-Y+. 154213393339382 0-16*FCTlJX) 

X=. 3610049480575200  02 

Y=Y+. 817182344342 072D-15*FCTICX) 

X=. 3234662915390470  02 

Y=Y*.3237801 65772 9270-1 3 *FCT1C XI 

X=. 2886210181632350  02 

Y=Y+.  9 799 3 7928 8 72 709 D-  12  *FCT  II X ) 

X=. 2562863602245920  02 

Y=Y+. 2305899491 891 34D-10*FCT1CX) 

X=. 2263  088901 31 9680  02 

Y-Y+. 428 138 29 7 1040930- 09*FCT1CX) 

X=. 1935586094033610  02 

Y=Y+ .635060222662581 D-08*FCTl C X 1 

X=. 1729245433671530  02 

Y=Y*. 76045678 7912 0780-07 *FC T l C X > 

X*. 1493113975552260  02 

Y=Y+.7416404578667550— 06+FCTIIX1 

X=. 1276369798674270  02 

Y=Y+. 5934541 61 286863D-05*FCT1 CXI 

X=. 10783018632 54 OOD  02 

Y=Y*. 392034196 798 7950-04 *FCT1 CXI 

X=. 8982940924212600  01 

Y=Y*. 21486491 8801 3640-03 *FCTl(X ) 

X=. 7358126733186240  01 

Y=Y+. 9808033066 14 9550-03 *FCT1 CX) 

X=. 5903958504174240  01 

Y=Y+. 373881 629461 1520-02*FCTICX) 

X=. 4616456769749770  01 

Y=Y+.l 19182 148348386D-01*FCTl CXI 

X=. 3492213273021990  01 


52 


AEDC-TR-74-51 


Y=Y+. 317609 12 5091 751D-01 *FCTl(X) 

X=. 2528336706425790  01 

Y=Y*. 7057  8623  8657 1740-01 *FCT 1 (X ) 

X=. 1722408776444650  01 

Y=Y+. 1299837862860720  00*FCT1(X) 

X=. 1072448753817820  01 

Y=Y+. 195903335972  88 ID  00*FCTltX) 

X=. 5768 846293 01 8 86D  00 

Y=Y«-. 2 35 2 1322966 9 34 80  00*FCT1(X) 

X». 23452610951961 9D  00 

Y=Y+.210443107933813D  00*FCT1U) 

X= .44439 365  3332 6 7 0D— 01 

Y=Y+. 1092183419523850  00*FCT1(X) 

RETURN 

END 
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TYPICAL  INPUT  DATA 


CARD  S 5 7 9 II  IS  IS  17  19  21  25  » 2T  29  SI  33  35  37  39  41  43  49  4T  49  SI 

# 

I 


2 


TYPICAL  OUTPUT  DATA 


4^  N«  5 14101)4-  0.66700  UF«  0.96800  DO-  0.10150  BO-  0.0 


'<40G4S0N*S  TRAJECTORY 


CllNST ANT  STRENGTH  JET 

BETA  XI 

0.0  0.0 

0.39589490  00  0.18585410-01 

0.10317460  01  0.148-8330  00 

0.13109660  01  0.50180600  00 

0.14223610  01  0.11894660  01 

0.14753850  01  0.2323176D  01 

P-0.200  0-0.400 


X 0W1  0W2  DU  1 0U2  DW  OU 

-0.20000000  01  -0. 65577770-01  0.82796130-01  -0.35621200-01  0.2051757D-01  0.17218360-01  -0.15103630-01 

-0.10000000  01  -0.18044750  00  0.14067690  00  -0. 15344770-01  0.33231800-01  -0.39770610-01  0.17887030-01 

0.0  0. 21100100  00  0.22049820  00  0.24154430  00  0.21752760-01  0.43149920  00  0.26329700  00 

0.10000000  01  0.90221030  00  0.25633430  00  3.25086670  00  - 0.51175370-01  0.1158545D  01  0.19969130  00 

0.20000000  01  0.12701150  01  0.21487160  00  -0.21800110  00  - 0.11388280  0 0 0.14849870  01  -0.33188380  00 


2ET4 

0.0 

-0.13340000  00 
-0.26630000  00 
-0.40020000  00 
-0.53360000  00 
-0. 66700000  00 
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NOMENCLATURE 


A 

Ai,  A2 

a 

Bi,  B2 
b 
C 
c„ 


-n 


'n 


n 


n 


*n 


En 

F 
F. 
f 

Gr 

h 

J.I.  P 
K 

Ke 

Ln 

-®n 


Mj 

m 

N 


n 


Area  of  jet 

Series  coefficients.  Appendix  B 
Width  of  a slot 

Series  coefficients.  Appendix  B 
Semiwidth  of  rectangular  wind  tunnel 
Cross-sectional  area  of  rectangular  tunnel,  4hb 
f cos  /3n 

Series  coefficient,  Appendix  B 
Slot  spacing 
q sin  /3n 

Initial  jet  diameter 
Series  coefficient.  Appendix  B 
Normalized  geometric  slot  parameter,  K/h 
Series  coefficient.  Appendix  B 

<q2+mV>1/2 

Series  coefficient.  Appendix  B 

Semiheight  of  rectangular  wind  tunnel 

Jet  intersection  point  with  tunnel  boundary 

Geometric  slot  parameter,  Eq.  (4) 

Equivalent  slot  parameter  for  quasi-linear  boundary 
condition,  2K 

Length  of  nth  straight  line  segment  approximation  to 
curved  wake 

Normalized  segment  length,  Ln/b 

Series  coefficient.  Appendix  B 

y- coordinate  Fourier  transform  parameter 

Number  of  straight  line  segments  approximating  the 
curved  wake 
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Nn 

P 


P 

Poo 

Pt 

L® 

Q 


q 

R 

R 


e 


R 


n 


r 


S 


U* 

uh 


X,  Y,  Z 
x,y,  z 

*0 

/3 


Series  coefficient.  Appendix  B 

Slot  parameter,  (1  + F)  P = 0 is  closed  wall  and 
P = 1 is  open  jet 

Series  coefficient.  Appendix  B 

Local  static  pressure 

Free-stream  static  pressure 

Free -stream  total  pressure 

Porosity  parameter,  (1  + l/R)-^,  Q = 0 is  closed  wall 
and  Q = 1 is  open  jet 

x- coordinate  Fourier  transform  parameter 
Porosity  parameter 

Pseudo  porosity  parameter  for  quasi-linear  boundary 
condition,  Eq.  (7) 

Series  coefficient.  Appendix  B 

Distance  from  any  field  point  to  element  of  vortex  ring 
representing  the  V/STOL  model 

Distance  along  the  jet 

Normalized  distance  along  the  jet,  S/b 

Initial  jet  velocity 

Free-stream  velocity 

Perturbation  velocity  in  axial  direction  induced  by 
horizontal  walls 

Interference  perturbation  velocity  in  axial  direction, 

a^i/dx 

Interference  perturbation  velocity  in  vertical  direction, 
positive  downward,  d(j>  j/9z 

Average  velocity  at  rotor  disk,  1/2  (dy/ds) 

Cartesian  coordinates  in  physical  dimensions 

Normalized  cartesian  coordinates,  X/b,  Y/b,  Z/b 

Zero-lift  angle  of  attack 

Inclination  angle  between  wake  trajectory  and  vertical  axis 
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&0  Initial  wake  angle 

r Circulation 

y Normalized  circulation,  T /U^b  ■ 

6U  Streamwise  interference  factor,  (Cui)/(Aw0) 

6W  Upwash  interference  factor,  (Cwi)/(Aw0) 

S Z location  of  jet  doublet 

d Vortex  angle 

A Tunnel  height -to -width  ratio,  h/b 

£ X location  of  jet  doublet 

<i>  Perturbation  velocity  potential 

0 Normalized  perturbation  velocity  potential,  (D/U^b 

0 k Normalized  perturbation  velocity  potential  due  to  the 

horizontal  walls 

(f>i  Normalized  interference  velocity  potential 

<£m  Normalized  perturbation  velocity  potential  due  to  the 

model 

0 v Normalized  velocity  potential  for  ah  infinite  row  of  model 

images 
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